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Abstract 


In this paper, we constrain the tilt of the power spectrum of relic gravitational waves by 
combining the data from BICEP2/Keck array and Planck (BKP) and the Laser Interferometer 
Gravitational-Waves Observatory (LIGO). From the data of BKP B-modes, the constraint on the 
tensor tilt is nt = 0 . 661(44 at the 68 % confidence level. By further adding the LIGO upper limit 
on the energy density of gravitational waves, the constraint becomes n* = — 0 . 761 q 52 at the 68 % 
confidence level. We conclude that there is no evidence for a blue-tilted power spectrum of relic 
gravitational waves and either sign of the index of tensor power spectrum is compatible with the 
data. 
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The primordial gravitational waves can be generated during inflation [iH^. The simplest 
inflation models show the power spectrum of primordial gravitational waves to be adiabatic, 
Gaussian and nearly scale-invariant. The primordial gravitational waves contribute to both 
the total intensity and polarizations for the anisotropy of the cosmic microwave background 
(CMB) j6-12|. Especially, the effects of primordial gravitational waves on the B-mode 
polarization might be detectable at the range I < 150. Recently, Background Imaging of 
Cosmic Extragalactic Polarization (BICEP2) 13|] reported an excess of B-mode power over 
the base lensed-ACDM expectation at the range 30 < £ < 150 with a significance of more 
than 5a. However, there were debates on whether the B-mode signals are resulted from the 
primordial gravitational waves or from the interstellar dust polarization 14l-ll7l|. Actually, 


the power of dust polarization has the same magnitude as the BICEP2 B-mode signals 
Most recently, a joint analysis of the data from BICEP2/Keck array and Planck (BKP) 


i 


19| 


yielded an upper limit on the tensor-to-scalar ratio, i.e., ro.os < 0.12 at the 95% confidence 
level (C.L.) which is compatible with the Planck temperature-only limit. The significance 
is too low to claim a detection of primordial gravitational waves. 

The tilt of power spectrum of relic gravitational waves is one other significant quantity. 
It describes the shape of tensor power spectrum. In general, the tensor tilt rit is defined by 


m) = A{- 

hjr 


( 1 ) 


where Pt{k) denotes the amplitude of the primordial tensor power spectrum at scale k, and 


kp = O.OlMpc”^ is the pivot scale. In the canonical single-field slow-roll inflation models, 
the tensor tilt is determined by a consistency relation between r and nt, namely, r = —8nt 
20|. Because of the upper limit ro.os < 0.12, the spectrum of relic gravitational waves 
is nearly scale-invariant. In the inflation model, the tensor tilt is generally predicted as 
Ut = —2e 20|, l2l|. The inflation requires a/a = H'^{\ — e) where e = H/H^, and thus 


—2 < Ui < 0. However, there are also certain alternatives of the inflation model, which 
predict different tensor tilt. For example, the ekpyrotic model j^| predicts a blue-tilted 
tensor power spectrum, i.e. n* = 2. Here we take the tensor tilt rit as a fully free parameter 
in our analysis. 

Besides constraining the tensor-to-scalar ratio, the current data can also constrain the tilt 
of the primordial tensor power spectrum. In this paper, we shall make a joint analysis of the 
CMB B-mode polarization and Laser Interferometer Gravitational wave Observatory (LIGO) 
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23| data to constrain the tensor tilt rit. The CMB B-niode polarization data coming from a 


joint analysis of the data BICEP2/Keck array and Planck 19| can be used to constrain the 
tilt of the power spectrum of relic gravitational waves at large angular scales. In this paper, 
we are just interested in the BB bandpowers 1 — 5 which are taken between BICEP2/Keck 
and the 217 and 353 GHz bands of Planck. As a complement, the data of LIGO place an 
upper limit on the energy density of relic gravitational waves at a specihc frequency band 
around lOOHz. 

The data of LIGO refers to the upper limit on the energy density of relic gravitational 
waves at the frequency band around lOOHz. The upper limit is given by Ggvk < 5.6 x 10“® 
at the 95% conhdence 


23| . In generah 


24 


26| 


wavenumber k is determined by 

^Gwik) = 


he energy density of relic gravitational waves at the 

Pt{k) 


12HI 


T{r]o,k) 


( 2 ) 


where rjo denotes the conformal time today, Hq is the Hubble parameter today, T{ri, k) is the 
transfer function of tensor perturbations, and the overdot denotes a cosmic time derivative 
d/dt. The wavenumber k of relic gravitational waves is related to the frequency / by / = 
k/^Ti. The tensor transfer function T{ri,k) describes the evolution of relic gravitational 


waves in the universe. It has an analytical approximation, namely, 

-3j2{kr]o)n, 
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T{r]o,k) = 


Vo 
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eq 
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( 3 ) 


where rjo = 1.41 x lO'^Mpc is the conformal time today, and keq = 0.073Gm,h^Mpc“^ denotes 
the wavenumber relating to the Hubble horizon at the time of matter-radiation equality. Here 
flm and h denote the matter density parameter and the reduced Hubble parameter today, 
respectively. Since LIGO is sensitive to the relic gravitational waves of the wavenumber 
k ^ keq, the energy density of relic gravitational waves today is given by [26 1 

15 VtLA.r ^ ^ 


^Gw{k) — 


r2^4 7. 


( 4 ) 


16 H^r]^klq \kp 

The data of LIGO can be a complement to the GMB data, since both observations refer 


to two very different cosmological sea! 
constraints on the inflationary physics 


es. Thus, one can make certain multi-wavelength 


We add a prior into the public Markov Ghain Monte Garlo sampler (GosmoMG) 29| to 
account for the LIGO upper limit on the energy density of relic gravitational waves. We 
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consider two combined datasets: one is the BKP B-mode data; the other one refers to a 
combination of the BKP data and the LIGO upper limit (BKP+LIGO). We just consider 
the BKP B-mode data with the bandpowers 1 — 5 (20 < £ < 200), which are insensitive to 
the lensing effects. Actually, the extra four bandpowers would lead little changes on our hnal 
results. The cosmological model considered here is the base lensed-ACDM model-|-tensor 
cosmology. The parameters of the base lensed-ACDM model are hxed as the “BKP cos¬ 
mology parameters” in the CosmoMC, namely, 1006 *mc 5 D lii(10^°As), n^) = 

(0.0220323,0.1203761,1.0411,0.0924518,3.16,0.9619123). Thus, the parameter space just 
includes the tensor-to-scalar ratio (r), the tensor tilt (rit). 

Our results are summarized in Table [H The marginalized contour plot and the likelihood 


parameter 

BKP 

BKP-hLIGO 

68 % C.L. 

95% C.L. 

68 % C.L. 

95% C.L. 

r 

< 0.055 

< 0.099 

< 0.059 

< 0.106 

nt 

0.66tj;«3 


''J-0.52 

-o.ietlf. 


TABLE I: The 68 % and 95% limits for the parameters r and nt from the BKP only and BKP-I-LIGO 
datasets. 

distributions of the parameters r and rit are illustrated in Fig. [H If the BKP dataset is 
considered, the constraints on the tensor-to-scalar ratio r and the tensor spectral index rit 
are r < 0.099 (95% C.L.) and rit = 0.66l};44 (68% C.L.), respectively. If the BKP-I-LIGO 
dataset is considered, the constraints become r < 0.106 (95% C.L.) and rit = 

(68% C.L.), respectively. Our results indicate that a scale-invariant power spectrum of relic 
gravitational waves is compatible with the data. Furthermore, the ekpyrotic model in which 
rit = 2 is ruled out at much more than 3a level once the LIGO dataset is taken into account. 
The canonical single-held slow-roll inhation is within la. Note that the errorbar of rit is still 
too large to test the inhationary consistency relation. 

In this paper, we make a constraint on the tilt of power spectrum of relic gravitational 
waves by combining the data of BKP and LIGO. The bounds on the tensor tilt are given 
by lit = 0.66^j;44 at 68% C.L. for the BKP dataset, and n* = —0.761 q; 52 at 68% C.L. for 
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FIG. 1; The marginalized contour plot and likelihood distributions of the parameters r and rit 
from the BKP only (red) and BKP+LIGO (blue) datasets. 


the BKP+LIGO dataset. We find that there is no evidence for the blue-tilted tensor power 
spectrum, and either sign of the tilt of tensor power spectrum is compatible with the data. 
The LIGO upper limit on the energy density of relic gravitational waves gives a stringent 
upper bound on the tensor tilt. However, the lower bound on rit is still quite loose. Actually 
the temperature auto-correlations are sensitive to the negative value of the tilt of tensor 
power spectrum 30, l3l| . Thus, the constraint on the negative part of the tensor tilt will be 


signihcantly improved once the the GMB TT spectrum is taken into account. Even though 
we obtain certain constraints on the tensor tilt just from GMB B-mode and LIGO datasets, 
it is worthwhile to further combining other datasets, such as the Planck TT spectrum, to 
make a more stringent constraint on the tilt of tensor power spectrum in the near future. 
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